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The wave function of the second + state of 12 C which was obtained long time ago by solving 
the microscopic 3a problem is shown to be almost completely equivalent to the wave function of 
the 3a condensed state which has been proposed recently by the present authors. This equivalence 
of the wave functions is shown to hold in two cases where different effective two-nucleon forces are 
adopted. This finding gives strong support for interpreting the second + state of 12 C which is the 
key state for the synthesis of 12 C in stars ( 'Hoyle' state ), and which is one of the typical mysterious 
+ states in light nuclei, as a gas-like structure of three a particles, Bose-condensed into an identical 
s-wave function. 



The a clustering nature of the nucleus 12 C has been 
studied by many authors using various approaches 0. 
Among these studies, solving the fully microscopic three- 
body problem of a clusters gives us the most impor- 
tant and reliable theoretical information of a clustering 
in 12 C within the assumption that no a cluster is dis- 
torted or broken except for the change of the size param- 
eter of the a cluster's internal wave function. As rep- 
resentatives for the solution of the microscopic 3a prob- 
lem where the antisymmetrization of nucleons is exactly 
treated, we here quote two works: one by Uegaki et al. 
and the other by Kamimura et al. both of which 
were published almost a quarter century ago. In these 
works, the 12 C levels are described by the wave function 
of the form A{x(s, t)(f>^} with A standing for the anti- 
symmetrizer, (f>^ = tp(ai)<p(a2)4>(ot3) for the product of 
the internal wave functions of 3 a clusters, and s and t 
for the Jacobi coordinates of the center-of-mass motion of 
3 a clusters. Here 4>((Xi) (i = 1, 2, 3) is the internal wave 
function of the a-cluster ctj having the form 4>{a.i) oc 
exp[-(l/8& 2 )Y^ m>n (r iTO - r in ) 2 ]. The wave function 
x(s, t) of the relative motion of 3 a clusters is obtained by 
solving the eigen-energy problem of the full three-body 
equation of motion; (<^|(i? - E)\A{x(s,t)(j)^}) = 0, 
where H is the microscopic Hamiltonian consisting of 
the kinetic energy, effective two-nucleon potential, and 
the Coulomb potential between protons. The difference 
between the works by Uegaki et al. and Kamimura et al. 
lies in the adopted effective two-nucleon forces, besides 
the differing techniques of solution. 

Both calculations by Uegaki et al. and Kamimura et 
al. reproduced reasonably well the observed binding en- 
ergy and r.m.s. radius of the ground 0^~ state which is 
the state with normal density, while they both predicted 
a very large r.m.s. radius for the second Oj" state which 
is larger than the r.m.s. radius of the ground OJ 1 " state 
by about 1 fm, i.e. by over 30%. The observed 0^ state 
lies slightly above the 3a breakup threshold and the en- 
ergies of the calculated 0^ state reproduced reasonably 
well the observed value although the value by Uegaki et 



al. is slightly higher than the 3a breakup threshold by 
about I MeV. The second + state of 12 C is well known 
as the key state for the synthesis of 12 C in stars ( Hoyle 
state ) and also as one of the typical mysterious + states 
in light nuclei which are very difficult to understand from 
the point of view of the shell model. For the understand- 
ing of the nature of the 0^" state with dilute density, the 
analysis by Uegaki et al. of the Reduced Width Ampli- 
tude (RWA) function of the 8 Be-a breakup is very useful. 
The RWA function which is defined as yh{p) — 

^12\/8W{M 8 Be,L)^a)Y L (p)}j T o\A{x(s,t)<Pl}) with 
p standing for the realtive coordinate between 8 Be and 
a, proved to have similar magnitude for all partial waves 
L ( L= 0, 2, 4) for the ground 0^ state but it turned out 
to be large only for L=0 for the Oj" state. This result for 
the 0^ state with dilute density implies that the 0^ state 
has a gas-like structure of 3 a-particles which interact 
weakly among one another, predominantly in relative S 
waves. This understanding of the 0^ state structure had 
been already presented by Horiuchi on the basis of the 3a 
OCM ( orthogonality condition model ) calculation Q, 
and is quite different from the picture of a 3a linear-chain 
structure [j| for this state. It should be mentioned here 
that both calculations by Uegaki et al. and Kamimura et 
al. reproduced well not only the energy but also other ob- 
served quantities related to the Oj" state indicating that 
their wave fuctions of the 0^" state are highly reliable. 
For example, the reduced a-decay widths of the Oj" state 
calculated by Uegaki et al. and Kamimura et al. at the 
channel radius a = 7 fm are 0.39 and 0.56, respectively, 
while the observed value is 0.38. The calculated values 
of the monopolc matrix element M(0j" — * 0^) by Ue- 
gaki et al. and Kamimura et al are 6.6 fm 2 and 6.7 fm 2 , 
respectively, while the observed value is 5.4 fm 2 . 

Recently, based on the investigations by Ropke, 
Schuck, and coauthors on the possibility of a-particle 
condensation in low-density nuclear matter [6], the 
present authors proposed a conjecture that near the na 
threshold in self-conjugate An nuclei there exist excited 
states of dilute density which are composed of a weekly 
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interacting gas of self-bound a particles and which can 
be considered as an na condensed state 0. This conjec- 
ture was backed by examining the structure of 12 C and 
16 O using a new a-cluster wave function of the a-cluster 
condensate type. The new a-cluster wave function actu- 
ally succeeded to place a level of dilute density ( about 
one third of ground state density ) in each system of 
12 C and 16 O in the vicinity of the 3 respectively 4 a 
breakup threshold, without using any adjustable param- 
eter. In the case of 12 C, this success of the new a-cluster 
wave function may seem rather natural because, as we 
explained above, we had already known that the micro- 
scopic 3a cluster models had predicted that the 0^ in 
the vicinity of the 3a breakup threshold has a gas-like 
structure of 3a-particles which interact weakly with each 
other predominantly in relative S waves. 

The new a-cluster wave function of the a-cluster con- 
densate type used in Ref. represents a condensation of 
a-clusters in a spherically symmetric state. The present 
authors extended the wave function so that it can de- 
scribe the a-cluster condensate with spatial deformation 
@- They applied this new wave function to 8 Be and 
succeeded to reproduce not only the binding energy of 
the ground state but also the energy of the excited 2+ 
state. In addition, they found that although the effect 
of the spatial deformation is not large, the introduc- 
tion of the spatial deformation brought forth a 100 % 
overlap of the condensate wave function with the "ex- 
act" wave function given by the microscopic 2a cluster 
model which solves the 2a-cluster equation of motion, 
(c/> 2 \(H-E)\A{x(r)(j) 2 a }) = 0. This fact forces us to mod- 
ify our understanding of the 8 Be structure from the 2a 
'dumb-bell' structure to the 2a dilute ( gas-like ) struc- 
ture. 

The purpose of this short note is to report on our study 
of 12 C using the extended 3a condensate wave function 
with spatial deformation and comparing the obtained re- 
sults for 12 C with those of the "exact" 3a cluster model 
wave functions by Uegaki et al. and by Kamimura et al. 
The most remarkable result of this comparison is that the 
O2 wave functions by Uegaki et al. and by Kamimura et 
al. are almost completely equivalent to our condensate 
wave functions with slight spatial deformation which are 
obtained by using the same effective two-nucleon force 
as Uegaki et al. and Kamimura et al., respectively. This 
result implies that the "exact" 3a cluster model wave 
functions for the second 0^" state of 12 C can definitely be 
interpreted as 3a-particle Bose-condensed state. 

The wave function of the na-cluster condensate with 
spatial deformation was introduced in Ref. [8| and the 
detailed explanation of it is given there. So here we give 
a brief explanation which is necessary in this paper. The 
wave function has the form 

$na[Px,0y,0z) = J d 3 Ri ■ ■ ■ d 3 R n 
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where Xi = (1/4) X} ra =i r in * s the center-of-mass coor- 
dinate of the zth a-cluster on, 4>(ai) is the same internal 
wave function of the a-cluster a^ as the previous micro- 
scopic 3a cluster model, B 2 = b 2 + 2fl 2 (k = x, y, z), 
and $ B (J?i, • • • , R n ) is Brink's a-cluster model wave 
function jj- It is to be noted that $ n a(Px, P v , fl z ) ex- 
presses the state where n a-clusters occupy the same spa- 
tially deformed center-of-mass orbit exp[— {2/B 2 )X 2 — 
{2/Bl)X 2 -(2/B 2 )X 2 }, while the internal a-cluster wave 
functions stay spherical. § n a{Px, Py , Pz) can be written 
as a product of the total center-of-mass wave function 
and the internal wave function $ n a(Px, fly, Pz) as 

V(A,A.A)«^{-^-^-^} 



X®na(Px,Py,Pz), 
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na (Px, 0y, 0z) = A [ exp { - £ ( & (Xix ~ X Gx f 

i=l x 

+ J2 (X iy -X Gy f + ^{X iz -X Gz f) }«/>(«!) • • • 4>(a n ) 



All the calculations are made not with & na (Px, P y , Pz) 
but with & na (Px, P y , Pz) which is an eigen state of to- 
tal momentum with eigen value zero. In this paper we 
assume axial symmetry of the deformation around the 
intrinsic z-axis and put fl x — fl y . The a-condensed wave 
function with good angular momentum which is obtained 
by spin projection is then written as 

®L(Px = Py,Pz)= fd cos 6d J M0 (9)Ry(e)$ na (P x =P y ,P z ), 



where R y {9) is the rotation operator around the intrinsic 
y axis which rotates & na by an angle 9, and d , [ I0 (9) is 
the small D-function. 

As effective two-nucleon forces, we adopt the ones of 
Uegaki et al. and Kamimura et al. One is the Volkov 
force No.l 01 with Majorana parameter M = 0.575, 
used by Uegaki et al., and the other is the Volkov force 
No. 2 [l(| with Majorana parameter M — 0.59, used by 
Kamimura et al. Hereafter, we refer to the former force 
as force I while the latter force is referred to as force II. 
We adopt the same values for the oscillator parameter b 
as Uegaki et al. and Kamimura et al., namely b = 1.41 
fm for force I while b — 1.35 fm for force II. 

In Fig.l we give the contour map of the J" = + bind- 
ing energy surface corresponding to the spin-projected 
state &ia°(flx = P v ,Pz) in the two parameter space, 
Px{= Py) and p z . The adopted effective force for this 
energy surface is force II. We see the energy minimum 
at P x (= P y ) — 1.5 fm and P z = 1.5 fm, which means 
that the minimum has a spherical shape. The minimum 
energy of —87.68 MeV is about 1.7 MeV higher than the 
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FIG. 1: Contour map of the energy surface of the + state 
in the two parameter space, (3 X (= @ y ) and (3 Z . The adopted 
effective force is force II. 

binding energy of —89.4 MeV obtained by Kamimura ct 
al. for the ground Of state. The energy surface in the 
case of force I is similar to the energy surface of Fig.l. 
The minimum energy obtained by the use of force I is 
-86.09 MeV and it is about 1.8 MeV higher than the 
binding energy of —87.92 MeV obtained by Uegaki et al. 
for the ground state. 

In Fig. 2 we give the contour map of the J v = + bind- 
ing energy surface corresponding to the state orthogonal- 
ized to the minimum energy state $3^°(min.) ( the state 
at the minimum energy point in Fig.l). The orthogo- 

nalized state is denoted as P±$ia = Py>0z) an d 1S 
expressed as 

= (1 - |$& J= °(niin.)> ($& J = ( min.)\)*£°(p x = P y ,p z ), 

$& J = (niin.) = $3 7 Q = (min.)/^/(^=0(min.)|%°(min.)) 

The adopted effective force for Fig. 2 is force II. We see 
an energy minimum at f3 x {= (3y) = 5.7 fm and (3 Z = 1.3 
fm in the oblate region of the map and a second energy 
minimum at (3 X (— (3 y ) — 2.9 fm and f3 z = 9.4 fm in the 
prolate region of the map. The minimum energy value 
is —81.55 MeV and, what is very remarkable, this value 
is almost the same as the binding energy of —81.66 MeV 
obtained by Kamimura et al. for the second 0^ state. 
The minimum energy of —81.55 MeV is close to the sec- 
ond minimum energy of —81.39 MeV, and there is a valley 
with an almost flat bottom connecting these two minima. 
An almost flat bottom of the valley means that the en- 
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FIG. 2: Contour map of the energy surface corresponding 
to the + state orthogonalized to the state at the minimum 
energy point in Fig. 1. The adopted effective force is force II. 



ergy of the spherical configuration is only slightly higher 
than that of the deformed configuration, namely the en- 
ergy gain due to the deformation is small. The energy 

surface by the orthogonalized state P±&ia°(Px = PyPz) 
in the case of the force I is similar to the energy surface 
of Fig. 2, and the minimum energy of the orthogonalized 
state is —79.83 MeV. Here again it is very remarkable 
that this value is almost the same as the binding energy 
of —79.3 MeV obtained by Uegaki et al. for the second 
state. 

The fact that for each case of the two different effective 
forces a single orthogonalized state Pj_$3= °(/3 x = P y ,P z ) 
yields almost the same energy as the "exact" energy of 
the O2 state obtained by solving a full three-body prob- 
lem of the microscopic 3a-cluster model, strongly sug- 
gests that the 0^ state wave function given by the mi- 
croscopic 3a-cluster model is similar to the rather simple 
state P±&ia °{Px = P y ,Pz) as long as the adopted effec- 
tive two-nucleon force is reasonably realistic. 

We also calculated the J 77 = 2 + energy surface corre- 
sponding to the spin-projected state &iu 2 (Px = P y ,P z ) 
for the forces I and II. The minimum energies for the 
forces I and II are obtained to be —83.61 MeV at f3 x {= 
f3 y ) = 1.30 fm and (3 Z = 0.35 fm, and to be -84.65 MeV 
at (3 X (— fly) — 1-50 fm and (3 Z = 0.35 fm, respectively. 
These minimum energy values for forces I and II are both 
higher by about 2 MeV than the lowest 2 + energies by 
Uegaki ct al. and by Kamimura et al., respectively, whose 
values are shown in Table I. 

We also performed the GCM ( generator coordinate 
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TABLE I: Comparison of the minimum energy of the spin-projected energy surface, GCM eigen energy, and the energy given 
by the full 3a calculation. Comparison is made for the Of, Of, and 2f states for two cases of the effective two-nucleon force. 
The energy surface for the Of state means that of the orthogonalized state P±Qia°{f3 x = /3 y ,f3 z ). Energies are in MeV. 
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method ) calculation for J T = + and 2 + by superposing 
®ia(Px = P y ,Pz) over various sets of {(3 X ,^ Z ); 

E $M = P' y ,P'MH - E k )\$( a {0 x = p y ,p z )) 

x//(/3„/3z) = 0. 

The adopted values of (3 X are (3 X = (i — 0.5) fm with 
i = 1 ~ 6, and those of (3 Z is (3 Z = (j — 0.5) fm with 
j = 1 ~ 8. Hence the total number of the adopted grid 
points (f3 x ,{3z) is 48. The calculated eigen energies of 
the Of, Of, and 2f states are given in Table I for the 
two forces I and II. We have checked the convergence of 
the calculation of the eigen energies by changing the sets 
of {(3 X ,(3 Z ) for the GCM calculation. We see in Table 
I that all the GCM eigen energies of the Of , 0^~, and 
2f states are almost the same as the energies of the mi- 
croscopic 3a-cluster model in both cases of forces I and 
II. Since the eigen energies obtained by solving the full 
three-body problem of the microscopic 3a-cluster model 
are the "exact" energies, we can say by using the mini- 
max theorem of the variational problem that this almost 
complete equivalence of our GCM energies with the exact 
energies means that our GCM wave functions of the Of, 
Of, and if states are almost equivalent respectively to 
the Of, Of, and 2+ wave functions of the microscopic 3- 
cv-cluster model in both cases of force I and II. In order to 
check further this almost complete equivalence, we give 
in Table II the comparison of the calculated r.m.s. radii 
and monopole matrix elements M (Of — > Of) between 
our GCM and the microscopic 3a-cluster model. We see 
nice agreement of the calculated quantities between our 
GCM and the microscopic 3cv-cluster model. In Table II, 
we see that the large r.m.s. radius of the 0^~ state is also 
predicted by our GCM as by the microscopic 3a-cluster 
model, but at the same time we see that the calculated 
value corresponding to our GCM is slightly larger than 
that of the 3a-cluster model. We think the reason is be- 
cause our wave function of the Of state which contains 
a large amount of the components of the 3a condensed 
wave functions &ia°(Px = P y ,P z ) with large (3 X and/or 
(3 Z so as to yield large r.m.s. radius may have a longer tail 
than the former 3a-cluster model. This possibly longer 
tail behavior of the GCM Of wave function may explain 



the slight underestimation of the monopole matrix ele- 
ment of the GCM versus the 3a-cluster model through 
the slightly enhanced mismatch between the Of and Of 
wave functions in the GCM case. 

The fact that the second Of wave function of the micro- 
scopic 3a-cluster model is almost completely equivalent 
to our GCM wave function of the second 0^~ state which 
has a very large r.m.s. radius or equivalently very dilute 
density is very important. Since our GCM wave func- 
tion of the Of state expresses the Bose-condensed state 
of 3a-clusters, as is clear from its large r.m.s. radius and 
from its functional form, we can say that the second Of 
wave function of the microscopic 3a-cluster model ob- 
tained long time ago underlines the fact that the second 
+ state of 12 C in the vicinity of the 3a breakup thresh- 
old has a gas-like structure of 3a clusters with "Bose- 
condensation" . 

Now we discuss the relation between our GCM wave 
function of the Of state which we denote as ^gcivKO^) 
and the orthogonalized state P±^ia°(Px — /3 y ,f3 z ) with 
minimum energy which we denote by ^±(0f). Although 
the energy of 1 &±(0f) is almost equivalent to that of 
'J'gcm(OJ) and also to that of the Of wave function of 
the microscopic 3a cluster model, we cannot simply con- 
clude that ^±(0f) is almost equivalent to ^gcm(O^). It 
is because ^±(0f) is not yet guaranteed to be orthogonal 
to the Of wave function. The orthogonality of 1 &±(0f) 
to $3^° (min.) which is the state at the minimum energy 
point of the energy surface is not the same as the or- 
thogonality to the Of wave function, and ^±(0^) may 
contain some amount of the Of wave function. We there- 
fore calculated the squared overlap value of the two wave 
functions, |(*±(0J)|*gcm(0^))| 2 . The obtained values 
are 0.95 and 0.97 for forces I and II, respectively. These 
large overlap values mean that the GCM Of wave func- 
tions are very similar to ^±(0f) in both cases of force I 
and II and hence verify our former statement that the Of 
wave function of the microscopic 3a cluster model is very 
similar to a simple state 1 &±{0f) so long as the adopted 
two- nucleon force reasonably describes the physics. 

We also studied the magnitude of the spherical con- 
densate component contained in our GCM Of wave func- 
tions. For this purpose, we first constructed the projec- 
tion operator P sp h onto the functional space S sp h spanned 
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TABLE II: Comparison of the r.m.s. radii i? rms and the monopole matrix element M(0^ — ► 0]*~) obtained by the GCM 
calculation with those by the full 3a calculation. Comparison is made for two cases of the effective two-nucleon force. R r ms are 
m fm, and M(0+ 0+) are in fm . 
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6.45 



2.40 
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by spherical condensate wave functions as P sp h = 
J2k l^ , s P h)( VI/ sphl' where \l>g ph arc the orthonormal basis 
functions of the space S sp h- ^s P h are constructed as fol- 
lows, 

Efe°(/?; = 0' y = P' z )\®i:°(f3x = y = P*))9 k m 

The calculated values of |(*GCM(0j)|P sp h|*GCM(0^))| 2 
are 0.92 and 0.91 for forces I and II, respectively. Of 
course, we checked the convergence of the calculation by 
changing the number of the adopted components \&J ph in 
P sp h- The large magnitudes of these values imply that 
'I'gcm(OJ) is mostly composed of the spherical conden- 
sate component by more than 91 %. At the same time we 
have to note that some amount (less than 9 %) of the de- 
formed component which is orthogonal to the spherical 



component is necessary in order to have quantitatively 
good reproduction of the observed quantities. 

We finally make a remark on the ground 0^~ wave func- 
tion. Since this state has a normal radius and density, 
three a clusters overlap strongly with each other in this 
state, which is totally different from the situation of the 
Oj state where the mutual overlap of three, or even two 
a clusters is small. Therefore even though the 0^~ state 
is well represented by a superposition of our condensate 
wave functions (we recall that our wave function contains 
the Slater determinant as a limit case), it does not mean 
at all that the state has an a condensation character 
which is only valid for the gas-like state of a-clusters. 

In summary, we have shown that the Oj wave func- 
tion of 12 C which was obtained long time ago by solving 
the full three-body problem of the microscopic 3a clus- 
ter model is almost completely equivalent to the wave 
function of the 3a condensed state. This equivalence has 
been shown to hold for two different effective two-nucleon 
forces. This result gives us strong support to our opin- 
ion that the 0^" state of 12 C has a gas- like structure of 
3a clusters with " Bose-condensation" . A more detailed 
report of the present problem will be given elsewhere. 
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